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ABSTRACT 



X 

■ We present N-body simulations of elliptical galaxy encounters into dry merg- 

ers to study the resulting unbound intergalactic stellar population, in particular 
that of the post-Main Sequence stars. The systems studied are pairs of spher- 
ical galaxies without dark halos. The stellar content of the model galaxies is 
distributed into mass-bins representing low- and intermediate-mass stars (0.85 - 
8 M0) according to Salpeter's initial mass function. Our models follow the dy- 
namical evolution of galaxy encounters colliding head-on from initial low-energy 
parabolic or high-energy mildly-hyperbolic orbits, and for a choice of initial-mass 
ratios. The merging models with initial parabolic orbits have M2/Mi=l and 
10, and they leave behind respectively 5.5% and 10% of the total initial mass 
as unbound stellar mass. The merging model with initial hyperbolic orbit has 
M2/Mi=l, and leaves behind 21% of its initial stellar mass as unbound mass, 
showing that the efficiency in producing intergalactic stars through a high-energy 
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hyperbolic encounter is about four times than through a parabohc encounter of 
the same initial mass ratio. By assuming that all progenitor galaxies as well as 
the merger remnants arc homologous systems we obtained that the intergalactic 
starlight is 17% and 28% of the total starlight respectively for the parabolic and 
hyperbolic encounters with M2/ Mi=l. In all models, different mass stars have 
the same probability of becoming unbound and feeding the intergalactic stellar 
population. 

Subject headings: Galaxies: elliptical and lenticular, cD; interactions; stellar 
content. Stars: AGB and post-AGB. Planetary nebulae: general. 

1. Introduction 

Intergalactic (IG) starlight, both diffuse and from resolved stars, has been intensively 
studied in the past decade in the IG medium of both groups and clusters of galaxies. The 
long-lived populous class of low- and intermediate-mass stars, in the form of red giants (RGB 
and AGB), planetary nebulae (PNe), and their diffuse hght, has been observed in the IG 
medium of different types of galaxy associations, from poor groups like the M81 group of 
galaxies (Feldmeier et al. 2004a) and the Leo region corresponding to the HI cloud (Castro- 
Rodriguez et al. 2003), to compact, tidal groups (White et al. 2003), up to the Virgo 
(Ferguson et al. 1998; DurreU et al. 2002; Feldmeier et al. 1998, 2003, 2004b) and Coma 
(Gerhard et al. 2005) clusters, as well as in higher redshift clusters (Zibetti et al. 2005). 

The observations are telling us that a considerable fraction of post-Main Sequence (PMS) 
stars in galaxy associations is found between galaxies. The contribution of the unbound stars 
to the total mass and light of the association varies greatly, depending mainly on the galaxy 
concentration of the considered cluster or group. Observations of poorly populated groups 
seem to indicate that the upper limit of the IG contribution to the total light is very low 
(up to 1.6%, Castro- Rodriguez et al. 2003), while surveys of IG PNe, AGB and RGB stars 
in the nearby clusters indicate that the intra-cluster (IC) starlight contribution is roughly 
5 to 20% of the total starlight (DurreU et al. 2002; Aguerri et al. 2005; Feldmeier et al. 
2004b), depending on the assumptions made on the stellar populations probed and on the 
completeness of the sample surveyed. The observed range of the fraction of IG starlight is 
also supported by the observations of diffuse starlight in clusters at z~0.25 (Zibetti et al. 
2005). 

A host of explanations for the origin of the IG starlight have been proposed in a variety of 
studies. The tidal interactions between galaxies have been explored in some depth by several 
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Authors, among whose Merritt (1983) and Moore et al. (1996). An alternative scenario that 
is well suited for the IC environment has been proposed by Muccione and Ciotti (2004), where 
the stellar stripping from galaxies is driven by interactions between the stellar orbits within 
the galaxies and the cluster tidal field. Among all possible explanations for the existence 
of the IG starlight, the scenario of the elliptical galaxy merging {dry merging) has been the 
least explored. 

In this paper we present numerical models of dry merging of galaxies with different 
initial mass ratios, aimed at describing possible scenarios for the production of the unbound 
stellar mass in the IG medium. Our models are crafted in particular to account for the PMS 
IG population that is produced by the merging of elliptical galaxy pairs. The red galaxy 
merging scenario may not be very common at the present time in the core of galaxy clusters, 
where the velocity dispersion are of the order of 1000 km s~^. Nonetheless, red mergers have 
very recently proved to be the common evolutionary path to field (van Dokkum 2005) and 
cluster (Tran et al. 2005) elliptical galaxies, and may be also the path to produce IG starlight 
in galaxy groups, or in clusters periphery, where the the velocity dispersions is typically much 
lower than in young galaxy clusters cores (Arnaboldi et al. 2004). Dry mergers have been 
recently studied by Gonzalez-Garcia & van Albada (2005ab), and have proved to preserve the 
properties of the elliptical galaxy fundamental plane (Gonzalez-Garcia & van Albada 2003; 
Nipoti et al. 2003), providing a further handle to study mergers and stellar populations. 

In Section 2 we present our models, including the type of numerical experiment per- 
formed, the astrophysical input parameters, the considerations on stellar populations, the 
initial conditions and constraints, and a description of the methodology and stability tests 
of the numerical models. Section 3 illustrates our results, with particular attention to the IG 
stellar population produced. The discussion is in §4, where we present a limited comparison 
with the observational data and the likelihood that the dry-merging scenario might account 
for the observed IG starlight. This paper represent a first attempt at modeling the IG pop- 
ulation with dry merging of galaxies, and the parameter space is thus exploratory. More 
detailed models, and a larger array of calculated observable parameters, will be presented in 
forthcoming papers. 
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2. Models 

2.1. Galaxies, Initial conditions 

We perform numerical models of three cases of dry mergers, with a choice of initial mass 
ratios M2/M1—I or 10 ^. The initial galaxy conditions are similar to those described by 
Gonzalez-Garcia & van Albada (2005a), where a systematic study of the encounters between 
two spherical systems without a dark matter halo was performed. These models without dark 
matter are a good first approximation to study the IG population in galaxy associations. In 
fact, Gonzalez-Garcia & van Albada (2005b) have shown that merging of elliptical galaxies 
with dark halos produces a luminous particle distribution that is very similar to that resulting 
from the merging of elliptical galaxies without dark matter (Gonzalez-Garcia & van Albada 
2005a). On the other hand, similar models of those presented here, but with the inclusion 
of dark matter will be performed in the future to confirm our results, since Gonzalez-Garcia 
& van Albada (2005ab) did not compare the unbound stellar population derived from the 
merging processes with and without the inclusion of dark matter. 

We use isotropic spherical Jaffe (1983) models as initial conditions for our experiments, 
and the algorithm developed by Smulders & Balcells (unpublished, see Gonzalez-Garcia 
& van Albada 2005a for a detailed description of the algorithm). The projected surface 
mass density (hereafter, surface density) of such models presents a slope that decreases 
roughly as i?^/^, which makes it a suitable representation for elliptical galaxies, although the 
central parts present a cusp. The distribution function presents an analytical solution for the 
collision-less Boltzmann equation, allowing the implementation of N-body initial conditions. 

In Table 1 we summarize the characteristics of the initial galaxy pairs. Column (1) gives 
the run identification code, where the number indicates the mass ratio, the lower-case letter 
denotes the impact parameter (h defines the head-on impact), and the capital letter indicates 
the energy of the orbit (Pfor parabolic, and ^for zero energy at infinity, or hyperbolic orbit). 
Hereafter, we will identify the encounter models with this code. Column (2) gives the initial 
mass ratio of the colliding galaxies, column (3) gives their initial separation, column (4) gives 
their initial relative velocity, and columns (5) and (6) give respectively the impact parameter, 
and the orbital energy of the initial setup. 

In Table 1, and later in describing the model parameters, we use the model units unless 
otherwise noted. We adopt non-dimensional units with with Newton's constant of gravity 
G=l. In each run, the theoretical half-mass radius of the Jaffe model, rj, and the total mass 



In this paper the suffix 1 always refers to the less massive galaxy of the merging pair. 
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of the less massive galaxy, Mi, are also set to 1. The models may be compared with real 
galaxies using the following scaling: 

[M] = Mj = 4 X 10^^ Mo, (1) 
[R] = rj = 10 Kpc, (2) 
[t] = 2.4 X 10^ yr. (3) 

By adopting these units, the velocity unit is: 

[u]=414kms~^ (4) 
Following Jaffe (1983) notation, the mass inside a radius r is defined as: 



M(r) = M (5) 

In run lOhP the galaxy model with mass M2 is a scaled up version of the model with 
mass Ml with ten times more particles, constructed following the scaling relation between 
mass and radius given by Fish (1964): 

where K is a constant. Following Jaffe's definition, and Fish's relation, the theoretical half- 
mass radius for the more massive galaxy is 3.162. 

In order to limit the calculation time, and yet preserving the physical significance of the 
results, we modify JafTc's models in such a way to obtain working galaxy models with finite 
radii. We impose a cut-off radius to all galaxy models, with Ri = 10 for the less massive 
galaxies. Such radius includes only 91% of the mass in the theoretical Jaffe's model, thus 
we need to re-scale the half-mass radius as to keep Mi = 1. The re-scaled half-mass radius 
for the less massive galaxy is then equal to 0.82. 

Model IhP is an equal mass encounter between two galaxies on a parabolic orbit, where 
Ml = M2 = 1. The centers of the two galaxies are placed at an initial distance of 4 Ri. 
Model lOhP is an encounter between two galaxies with a mass ratio 10 on a radial parabolic 
encounter, initially placed at a distance 3 Ri + R2- Ri and R2 are the galaxy radii of the 
working models. Model IhZ is an encounter between two galaxies with Mi = M2 = 1, 
placed initially on a mildly-hyperbolic orbit. 
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The choices of initial conditions for the galaxies are adequate to represent the observed 
mergers that occur in elliptical galaxies. The galaxy mass ratios chosen for the simulations, 
M2/Mi=l and 10, are the extremes of the observed merging galaxy sample by van Dokkum 
(2005). 



2.2. Stellar components 

We populate our model galaxies with model stars that represent stellar masses in the 
m = 0.85 — 8 Mq mass range ^, whose progeny includes RGB, AGB, and PNe. Since 
the main scope of our modeling is to study the stellar population that produce PMS stars, 
we neglect all stars outside this mass range. In order to model these stars we use test 
particles whose masses are proportional to the stellar masses they represent. We assume the 
validity of Salpeter's (1955) initial mass function (IMF) \I'(m) oc m^^'^'^\ and consider three 
representative mass bins respectively for the progenitors of massive (3-8 Mq), intermediate 
(1.4-3 M0), and low-mass (0.85-1.4 Mq) PMS stars. The population of each mass bin 
corresponds to the integration of the Salpeter's mass function in that bin, scaled to the 
entire population considered, and ignoring stars outside the mass range. We then calculate 
the mass fraction for each bin as: 

/•mmax 

^bin = / m *(m) dm, (7) 
where mmin and mmax are the limits of the mass bin considered, and 

rS.O 

A= ^{m)dm^ 0.8777 (8) 

io.85 

is the normalization of Salpeter's law in the entire mass range. 

To characterize the mass bins in the N-body simulation we use a representative mass 
for each bin (1, 2, and 6 Mq). The model stars in the model galaxy with M = 1 have such 
masses that, after accounting for the Salpeter's IMF, the total galaxy mass is equal to unity. 

In Table 2 we summarize the characteristics of the stellar population in the M — 1 
galaxies. Columns (1) gives the mass bin in solar masses, column (2) gives the representative 
mass in that bin, colum (3) gives the mass fraction in that mass bin, column (4) gives the 
number of particles in the bin, and column (5) the (dimensionless) stellar mass for each 



^Hereafter, m refers to the stellar mass, to distinguish it from M, the galaxy mass 
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particle in that bin. Note that the mass of each particle is different from the mass of the 
star it represents, but the ratio between particle masses in different bins is the same than 
the ratio of the representative masses in that bin. This is a good representation, since we 
are interested in relative results for the mass bins. In runs IhP and IhZ both galaxies have 
number of particles per mass bin as in Table 2. The least massive galaxy of run lOhP also 
has stellar population as in table 2, while the more massive galaxy has an initial set up with 
ten times more particles in each mass bin, but with the same mass per particle. 

The initial distribution of the populations of particles with different mass is such that 
each stellar population follows the same Jaffe law. In this way, there are fewer particles 
from the high mass bin at any given radius. We use the distribution function (DF) as in 
Gonzalez-Garcia & van Albada (2005a), extending it for each mass bin. 

2.3. Integration method and stability 

We have used the parallel tree-code GADGET-1 (Springel et al. 2001) on the Beowulf 
Gluster at the lAC, where a typical run on 16 CPUs takes of the order of 1.5 x 10^ seconds. 
Gravitational Plummer Softening (sec e.g. Binney & Tremaine 1987, Eq. 2-194) was set to 
~ 1/10 of the half-mass radius of the less massive galaxy, with softening parameter e — 0.075. 
The tolerance parameter (see Barnes & Hut, 1986) was set to ^ = 0.8. Quadrupole terms 
were included in the force calculation. GADGET-1 uses a variable time step (ig), and we 
choose ts oc l/|a|°'^, where a is the particle acceleration. Wc set the maximum time step to 
1/100 of the half- mass crossing time, and the minimum time step to zero. 

We have checked the stability of our input initial model for 33 crossing times. In Figure 
1 we show the results from a test run that has been used to check the stability of our models. 
Therein, the evolution of the mass inside different radii has been followed at different mass- 
fraction levels, showing stability in a large range of t/ra-, where r^r is the half-mass crossing 
time. The test shows that the system relaxes for about 4 crossing times and remains stable 
thereafter. This initial relaxation is probably due to the presence of the particle softening 
in the code. 

Models were evolved for at least 10 dynamical crossing times of the merged system after 
merging, to allow the system to relax (reach virialization). Conservation of energy is sound 
in all the runs, with variations lower that 0.5% of the total energy. 

We performed an additional run of a model identical to IhP, except with twice as many 
particles in each mass bin. Our merging simulation on this additional model results in an 
intergalactic mass fraction of 5.22%, 5.24%, and 5.34% respectively for particles in the first. 
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second, and third mass bin respectively, proving that our results are stable against mass 
resolution. 



3. Results 

The three models of dry merging presented in this paper have different initial conditions 
both in the mass content of the parent galaxies and in their relative velocities, as described 
in §2. As model IhP is let evolve, the system passes through pericenter for the first time at 
t=85, when an exchange of orbital energy into internal energy occurs. The particles having 
initial binding energy close to zero gain enough energy at this stage to be expelled from the 
system and become unbound particles. The two galaxies finally met again 30 time units 
after the first encounter. This is the time when the actual merging occurs. New particles 
are able to become unbound at this time, as a consequence of new particle encounters and 
exchange of energy. The merging time is 2.8 x 10^ yr after the initial placement in orbit, 
with the unit convention given in §2.1. A very similar situation occurs in model lOhP, where 
the two unequal galaxies meet at pericenter after t=70 since the initial placement in orbit, 
then meet again 131 time units after the first encounter, with merging time of 4.82 x 10^ 
yr. The galaxies of model IhZ, initially in mildly-hyperbolic orbits, meet for the first time 
at t=34.5 and after 65.5 time units the second and final encounter occurs, ending in the 
merging episode. The merging time for IhZ is then 2.4 x 10^ yr. All runs are stopped after 
ten merging crossing times, after virialization has been reached. 

In Table 3 we give the characteristics of the stellar populations of our models after the 
merging has occurred. We give the run code (column 1), the mass bin (column 2), and the 
fraction of resulting unbound mass (column 3). We also give the unbound starlight fractions 
in columns (4) and (5), as described in the Discussion. 

By inspecting the data of Table 3 we find several interesting results from our merging 
simulations. First, we sec that for all runs the percentage of particles that become unbound 
does not depend on the mass bin. This result is expected, since the particle expulsion is due 
to gravitational acceleration, and indicates that if we were able to observe all IG stars with 
the same probability, we should recover the IMF of the original galaxies. Second, the final 
fraction of unbound mass is different for the three merging models: in IhP the fraction of 
unbound to total mass is Mjc/Mtot = (5.48 ± 0.1) x 10~^, which means that ~ 95% of the 
total initial mass still remains in the merger (Mm = 1.89 with Mtot = 2). In run lOhP the 
fraction of unbound to total mass is almost twice as high, with Mic/Mtot = (9.97±0.4) x 10~^, 
corresponding to a merger with 90% of the initial mass (Mm = 9.9 for Mtot = H)- Naturally, 
the larger fraction of unbound mass in model lOhP compared to IhP could be imputed both 
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Fig. 1. — A test run was made to check the stabihty of our initial models. The plot shows 
the evolution of the mass inside different radii, the top line gives the 99% mass radius, while 
the third hue from the bottom gives the 10% mass radius. 
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Fig. 2. — Binding energy versus distance (from the merger's center) of the model stars in 
the IhP model. The panels, left to right, represent the situation in the low, intermediate, 
and high-mass bins. 
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Fig. 3. — Same as in Figure 2, but for the lOhP model. 




Fig. 4. — Same as in Figure 2, but for the IhZ model. 
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to the initial mass ratio, but also to the high kinetic energy due to the higher initial relative 
velocity in model lOhP. The third model, IhZ, has the highest efficiency in producing 
unbound stars. The IG stellar mass in this case represent 21% of the total initial mass, and 
the merger retains only 79% of the initial mass. In this case, the merger has — 1.58 
from an initial total mass Mtot = 2. 

In Figures 2, 3, and 4 we show the location of the model stars in the radius ~ binding 
energy plane for each of the mass bins, for the 1/iP, lOhP, and IhZ models respectively. 
Both R and -Ebind in the plot are dimensionless variables, R being the distance from the 
merger's center. Note that the plots do not include the complete range where particles are 
found, rather the range where most particles are found. For example, out to a radius R=50 
we find 90% of the IhZ particles, and up to 99% of the IhP particles. These figures give 
a clear idea on the nature of the merger, and on the distribution of bound and unbound 
particles for each mass bin in each model considered. Unbound particles are found out to a 
distance of Rmax ~160, 600, and 200 respectively for models IhP, IQhP, and IhZ. We do 
find a mild anticorrelation between the particle mass and Rmax- By running a model similar 
to IhP but with twice as many particles per mass bin, we proved that the anticorrelation is 
due to the relative populations of the bins, and it does not have physical significance. 

In Figure 5, 6, and 7 we plot the surface density calculated within annuli of increasingly 
larger radii from the merger's center in models IhP, lOhP, and IhZ respectively. In these 
figures we plot the logarithm of the surface density fi against R^^^, where R is the distance 
from the merger's center. The three lines represent the different mass bins of stellar popu- 
lations. Note that these plots are non cumulative, and that the surface density have been 
evaluated in concentric annuli equally spaced in logR. This representation is useful to show 
where the galaxy merger profile dominates, and where the de Vaucouleurs (1959) profiles 
start to be perturbed by the unbound particles. By comparing the IhP (Fig. 5) and IhZ 
(Fig. 7) models, we see that the surface density profile is similarly perturbed, but the IG 
component outside the R^^^=2 annulus is much more important in the IhZ than in the IhP 
model. The de Vaucouleurs' slope is affected by the IG particles outmost of R^^'^ ~ 2.5 in 
the lOhP model (Fig. 6). In all models, test particles from the different mass bins contribute 
to the merger and the IG mass in very similar fashion. 



4. Discussion 

Elliptical galaxy encounters were modeled with N-body simulation, and they produce 
mergers in a fraction of the Hubble time. During the process they expel a fraction of their 
stellar content. The model stars expelled continue to evolve in the IG medium, eventually 
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Fig. 5. — Surface density versus i?^/^ in model IhP. The three hues represent the mass bins, 
where the sohd hne represents the low mass bin, the dashed line the intermediate mass bin, 
and the dash-dot line the high mass bin. 
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Fig. 6. — Same as in Figure 5, for model lOhP. 
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Fig. 7. — Same as in Figure 5, for model IhZ. 
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going through the red giant, AGB, and post-AGB phases. The fraction of unbound mass 
resulting through this process depends on the properties of the encounter. By comparing the 
two models with the initial parabolic orbits and different mass ratios we conclude that the 
model with the higher mass ratio produce a larger fraction of unbound mass with respect to 
the total (as well as the merger's) mass. By comparing the M2/M1 = 1 runs we see that the 
initial hyperbolic orbit results in a unbound mass fraction that is almost four times that of 
the parabolic orbit encounter. 

Since the elliptical mergers belong to the fundamental plane of elliptical galaxies (Gonzalez- 
Garcia & van Albada 2003), in the same way than the two original merging galaxies, we can 
determine the fraction of unbound-to-total starlight for all our models, simply by assum- 
ing that both the initial galaxies and the final mergers are homologous galaxies. If we 
assume that every galaxy of the merging pairs and all mergers obey the mass-luminosity 
relation M/L oc M°-^ (Jorgensen et al. 1996) we can calculate Lio/Liot, the ratio of the 
IG starlight over total starlight. In dry merging encounters the relation linking the total to 
individual galaxy luminosity can be written as Liq = {Li + L2 — L^). By using the mass 
fractions of unbound mass from Table 3 we found [Liq/ Ltot]ihP='^7%, [-^^iG/-^tot]iofeP=14%, 
and [LiG/i^tot]i?iZ=28%. The fractions of the unbound to merger's starlight are respectively 
[Lig/LJi/ip = 20%, [LiG/L^]whP = 17%, and [LiG/L^]ihz = 38% (see also Table 3). 

By recalling that our model galaxies are populated by low- and intermediate-mass stars, 
their Liq/L^ fractions may be compared with the observed PMS populations. The IG 
starlight ratios that we recover in our simulations are in broad agreement with the ob- 
servations of several IC and intra-group stellar populations. Following, we examine a few 
particular cases. 

The most studied case of IG stellar population is that of PNe in Virgo. The mechanism 
proposed in this paper is probably more likely to occur as a possible origin of the IC stellar 
population, or parts thereof, in the cluster periphery rather than at the center of a cluster 
such as Virgo. Arnaboldi et al. (2004) have studied the velocity dispersion of a sample of 
Virgo IC PNe and found that several fields have dispersion velocities much lower (247 km 
s~^) than the canonical Virgo dispersion of 800 [km s~^], obviously a consequence of the fact 
that the cluster is young and highly nonuniform. Feldmeier et al. (2004b) found that 16% 
of the starlight in the Virgo cluster is in the IG medium, independent on the location within 
Virgo, and this number is encompassed by our results. 

Nonetheless, the observed ratio of unbound to total starlight in the Virgo cluster derived 
from PNe must be used cum grano salis. To derive this ratio from observations one must 
evaluate the theoretical luminosity-specific PN density, «pn = B tpN (Renzini & Buzzoni 
1986) , which is based in the fuel consumption theorem for PMS stars. The question is whether 
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this theorem is adequate to describe a nebular population, where tp^, the PN lifetime, is 
not a stellar evolutionary time, determined by fuel consumption, but rather a time-scale 
depending on the hydrodynamic evolution and the photo-ionization on the nebulae. While 
tpN = 25, 000 yr is typically adopted, the correct lifetime for PNe to be observable at high 
luminosity is probably much lower, as the hydrodynamic models by Villaver & Stanghellini 
(2005) have shown. Furthermore, tpN has not been parametrized for the mass and chemistry 
of the progenitor stars. Feldmeier et al. (2004b) advise to use ccstars from other stellar 
sources, such as Virgo red giants (Durrell et al. 2002). 

The red giant population of the Virgo IC, first observed by Ferguson et al. (1998), 
accounts for approximately 10% of the cluster (evolved) stellar mass. Ferguson et al. (1998) 
indicated that the IG population is likely to originate from eUiptical and SO galaxies, for their 
higher frequency in the cluster and their older stellar populations. Intergalactic starlight has 
been observed in Fornax both in the form of PNe and other stars. Given the difficulty of 
PNe lifetime scaling, we prefer to use the results from IG novae observations (Neill et al. 
2005), indicating that ~16 to 41 % of the starlight in the Fornax cluster comes from IG 
stars. Both the Virgo and Fornax PMS IG stars might have dry merging origin, and these 
percentages are clearly in the range of our results (see Table 3, column 4). Group IG 
populations are typically evaluated to be 10 times lower than their cluster counterparts, 
although spectroscopic confirmation of, for example, the M81 PNe has not yet been published 
(Feldmeier et al. 2004a) and final counts are not available. 

The mere existence of low-mass PMS stars in a given stellar population requires a very 
old stellar population. Following Maraston's (1998) prescription, the turnoff mass in the 
first star mass bin correspond to ages in the ~4.8-22 Gyr range. If, for sake of reasoning, 
we assume that the progenitor galaxies were just formed at the time they were put in their 
relative orbits, none of their low-mass (0.85-1.4 Mq) stars would have reached the turnoff 
by the time the merger is completely formed. More realistically, the merging galaxy pairs 
contain aging stellar populations when their first encounter occurs. Accordingly to van 
Dokkum and collaborators (van Dokkum et al. 1999, van Dokkum 2005), dry merging could 
have been important at intermediate redshifts when clusters were still assembling. Therefore, 
part of the IC stellar population may have been expelled at those cosmic times. 

Zibetti et al. (2005) found that the intracluster light at large cluster radii is largely 
dominated by surface brightness excess around galaxies. From our Figure 5 we find that the 
IC contribution to the profile is above the R^^^ law from i?^/^ 1.5 to 2. This implies that 
the IC stars contribute the most to the surface density out to ftilOO Kpc from the merger edge 
of model IhP (or, with similar reasoning, about 60 Kpc from the merger edge in model IhZ) 
with the usual unit conversion, consistent with what observed by Zibetti et al. Naturally 
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this analysis is sensitive to the choice of units that we use to compare the models with real 
galaxies, as described in §2, but it is worth showing that our equal-mass encounter models 
are at least broadly consistent with the observations. 

In this paper we have compared our models with the observed IG populations. In 
other words, we have implicitly assumed that all IG stars derive from dry merging, and 
that all galaxies in a given observed cluster or group has gone through at least one merging 
process. In this extreme assumption the models correctly predict the observed IG starlight. 
Naturally, the case might well be that dry merging occurs only in a fraction of cases, and 
that other mechanisms arc at work in explaining the observed unbound starlight in galaxy 
associations. While the dry merging scenario is certainly helpful to account for a fraction 
of the observed IG light, it may not work well near cluster centers and in high velocity 
environments, where other mechanisms such as mass stripping due to hyperbolic encounters, 
or disk galaxy merging, might be more efficient. 

In the future, we plan to model similar encounters as those presented here, but including 
dark halos in elliptical galaxies. Furthermore, other phenomena related to dry merging, and 
a larger range of initial conditions, will be considered in future studies. 

We warmly thank Dr. Luca Ciotti for important comments on the manuscript and for 
scientific discussions, and Drs. Mark Dickinson and John Feldmeier for scientific discussion 
and bibhographic indications. L. S. thanks the European Southern Observatory and the I AC 
for their hospitality during the summers of 2004 and 2005, when this work was conceived 
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Table 1. Model input, galaxy pairs parameters 



Run M2/M1 n Vi b Eorb 



IhP 1 40 0.316 

10/iP 10 61.62 0.604 
IhZ 1 40 1.048 0.250 
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Table 2. Model input, stellar parameters 



bin [M0] 


M[M0] 




^particles 


^particle 


0.85 - 1.4 


1 


0.5151 


156060 


3.3 X 10-^ 


1.4 - 3.0 


2 


0.3443 


52121 


6.6 X IQ-^ 


3.0 - 8.0 


6 


0.1406 


7050 


2.0 X 10-^ 


total 




1 


215231 


1 
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Table 3. Model results 



Run 


bin 


MiG 


Lig/ Ltot 


Lig/ Lra 




[Mo] 


% 


% 


% 


IhP 


0.85 - 1.4 


5.53 






IhP 


1.4-3.0 


5.55 






IhP 


3.0 - 8.0 


5.37 






IhP 


total 


5.48 


17 


20 


lOhP 


0.85 - 1.4 


9.87 






lOhP 


1.4-3.0 


10.20 






lOhP 


3.0 - 8.0 


9.57 






lOhP 


total 


9.97 


14 


17 


IhZ 


0.85 - 1.4 


20.6 






IhZ 


1.4 - 3.0 


20.6 






IhZ 


3.0-8.0 


20.9 






IhZ 


total 


20.7 


28 


38 



